The Role of Thermal Phase Fluctuations in Underdoped YBa2Cu307_5 Films 
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The effect of thermal phase fluctuations (TPFs) on the afc-plane penetration depth, A(r), of 
thin YBa2Cu307_i (YBCO) films is found to be much smaller than expected from the paradigm 
of cuprates as weakly-coupled 2D superconducting layers. A 2D vortex-pair-unbinding transition 
is observed, but the effective thickness for fluctuations is the film thickness, not a CuO bilayer 
thickness. In a strongly underdoped YBCO film, Tc = 347^, TPFs suppress Tc by only about 3 K. 
They cannot be a significant factor in the suppression of Tc and emergence of the pseudogap with 
underdoping. 

PACS numbers: PACS numbers: 74.25. Fy, 74.40.-|-k, 74.76.Bz, 74.72. Bk 



Understanding the origin of the pseudogap in under- 
doped cuprates is one of the central issues in cuprate 
superconductivity. One proposal is that the tempera- 
ture, T* , where the pseudogap appears is actually the 
mean- field superconducting transition temperature, Tco, 
and some sort of superconducting fluctuation suppresses 
the onset of phase coherence down to the measured Tc, 
which vanishes at strong underdoping. It has been pro- 
posed that thermal fluctuations in the phase of the su- 
perconducting order parameter are an important effect 
in cuprates P, 0|- There are two good reasons to 
expect strong thermal phase fluctuation (TPF) effects 
in underdoped cuprates: first, coupling between Cu02 
(bi) layers is so weak that fluctuations should be quasi- 
two-dimensional; and second, the superfluid density ^5 
is small. These factors grow with underdoping. 

It is important to know whether TPFs alone can ac- 
count for a significant suppression of Tc- While much of 
our understanding of TPFs comes from studies of arrays 
of Josephson-coupled grains, 0j 13 ^ significant effort 
has gone into more realistic models 0,13 ■ Curty and 
Beck 6] find that fluctuations in the amplitude of the 
order parameter are more important than phase fluctu- 
ations. In the present paper, we look for experimental 
evidence that thermal phase fluctuations, as understood 
from simulations of arrays of Josephson-coupled grains, 
play a significant role in suppressing Tc of heavily un- 
derdoped cuprates. 

The feature that we identify as the onset of strong 
TPFs is a rapid increase in downward curvature in ns(T), 
as found in numerical simulations of 2D and quasi- 
2D 1^ superconductors. These simulations find that 
when TPFs are strong enough to suppress ng about 
35% below its mean-field value, nonlinear effects (vortex- 
antivortex pairs; vortex loops) come into play, and as T 
increases further these nonlinearities rapidly suppress ns 
to zero. Not surprisingly, nonlinear effects emerge at the 
vortex-pair-unbinding transition temperature, Tlj*^, of a 



single layer. The downward curvature of ns{T) begins 
to grow at T w T*jj and diverges as (Tc - T)~^/^ at 
a 3D-XY transition 0, 0|. We are interested only in 
finding the onset of strong fluctuation effects, not details 
of the critical region, which are often obscured by inho- 
mogeneities anyway. The abrupt downturn in ns should 
be enhanced in real samples relative to quasi-2D classical 
calculations because of two effects: first, as discussed 
below, TPFs are suppressed for T below an effective De- 
bye temperature, so when they turn on, they turn on 
more rapidly than found in simulations; and second, to 
map simulations onto data, the Josephson coupling en- 
ergy in the simulations must be proportional to the mean- 
field superfluid density and hence should decrease as T 
increases. 

Quantitatively, the vortex-pair unbinding tempera- 
ture, T2D, of a 2D superconductor, thickness d, is pre- 
dicted by the well-known relation 18] : 
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where = d\^^ oc nsd is the 2D magnetic penetra- 
tion depth, and $0 is the flux quantum. Note that the 
superfluid density is proportional to AJ^. The hypothet- 
ical 2D transition temperature, 7^2*0 ' f'^'' ^ single Cu02 
bilayer in YBCO follows from Eq. (QJ with d = 1.17 nm, 
the center-to-center spacing between bilayers. In effect, 
Eq. predicts a transition when the thermal energy, 
ksT, equals the superconducting condensation energy in 
a characteristic volume 27r^^d. Measurements of ns in 
2D films of conventional s-wave superconductors validate 
this equation (see, e.g., ref. jl^). 

To locate the onset of TPF effects, we fit a quadratic 
(constant curvature) to ns{T) just below T^jj, then find 
the temperature where ns(T) drops below the fit. We 
first show that in a very thin (4 unit cells) optimally- 
doped YBCO film a rapid downturn in ns is clearly 
present, so there is nothing intrinsic to YBCO that pre- 
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eludes the usual vortex-pair unbinding transition. In the 
4-unit-cell-thick film, the observed Tc of 70 K is seen to 
be about 15 K below its mean-field value, Tqo « 85K. 
Thicker (8 and 10 unit cells) optimally-doped YBCO 
films also show downturns in ns, but at temperatures 
that suggest that the effective thickness for TPFs is 
the film thickness. A severely underdoped YBCO film, 
Tc = 34 ± 2K, also has a rapid downturn in ns at a 
temperature well above T^jj. The data show that Tc is 
suppressed below Tco by only 3 K. 

Our results seem to conflict with experimental support 
for strong TPF effects in YBCO provided by observa- 
tions of a critica l reg ion about 5 K wide in the superfluid 
density 0, 0, JSjand other properties of very clean 
YBCO crystals [lj,|l3- We argue that there is something 
anomalous about these results, however, since they are 
apparently very sensitive to disorder, and are thus not 
supported by data on slightly less clean YBCO, on YBCO 
films 19], or on BSCCO (see below) that is much more 
anisotropic than YBCO. 

This study reports superfluid density measurements on 
four samples: films A, B, and C were grown epitaxially 
by pulsed laser deposition (PLD) on NdGaOa substrates, 
as detailed elsewhere j^^l- They are fully oxygenated. 
Each 1 cm x 1 cm film consists of buffer layers of semi- 
conducting Pro.6Yo.4Ba2Cu307_5 that are 12 unit cells 
thick above, and 8 unit cells thick below, the YBCO film. 
The underlayer lessens the strain of substrate lattice mis- 
match on the YBCO, while the capping layer protects the 
YBCO film from damage during handling. Films A, B, 
and C are nominally 4, 8, and 10 unit cells thick; however, 
the top and bottom YBCO layers may not be perfectly 
smooth or homogeneous. Our conclusions are insensi- 
tive to this uncertainty. Film D was grown by PLD on 
a SrTiOa substrate at 760 C, then annealed at 600 C 
in a low pressure of oxygen so that it would be severely 
underdoped. There were buffer layers of PrBa2Cu307_5 
above and below the 40-unit-cell thick YBCO film. 

We measured A]^^(r) ~ dfiim/-^^ with a two-coil 
mutual inductance technique described in detail else- 
where ll^. Each film was centered between two coils 
roughly 2 mm in diameter and 1 mm long. In a typi- 
cal measurement, the sample was cooled to 4.2 K, and a 
current of roughly 100 was driven at 50 kHz through 
the coil pressed against the back of the substrate. In 
this geometry, the induced ac electric field in the film 
was parallel to the plane of the film and had azimuthal 
symmetry. It was very nearly uniform through the film 
thickness because film thicknesses were much less than 
the magnetic penetration depth, A. Thus, the conduc- 
tivities of all layers in the film were in parallel, and the 
measurement yielded the sheet conductivity of the entire 

film: CT(w,T)dfiim = tTi(w,T)£ifilm - jo-2(w,T)£ifiim- As T 

slowly increased, the voltage induced in the secondary 
coil, which was pressed against the film, was measured 
continuously. A]|^(T) was obtained with an accuracy of 




40 60 
T[K] 



FIG. 1: A^^ = dfiim/A^ vs. T measured at 50 kHz for 
optimally-doped YBCO films A, B, and C (thick solid lines). 
Vertical dotted lines locate the peaks in ai(50 kHz,T), which 
very nearly coincide with T2D (open circle). Thin curves are 
quadratic fits to A^^(r) just below T2]j (open circle). Insets 
enlarge the transition regions. 



about 3% from cr2rffiim by using the relation: 

Xj^iT) = ^0WCT2(T)dfilm, (2) 

where fxo is the permeability of vacuum. Uncertainty in 
rffiim enters only in calculating the 3D penetration depth, 
A~^(r), or the 2D penetration depth of a single unit-cell 
layer. 

Figure 1 shows X^^{T) vs. T (thick curves) measured 
for the thin, optimally- doped YBCO films. A2^(T) is 
quadratic in T at low T presumably due to small disor- 
der and a d-wave superconducting gap. The value of 
A]]^(0) for the 10 unit-cell-thick film is what we rou- 
tinely observe in thick YBCO films made by pulsed 
laser deposition. A]]^(0) decreases for thicknesses less 
than 10 unit cells, for reasons yet to be determined. 
Tc, denoted by vertical dotted lines, is defined as the 
center of the fluctuation peak in cri(50 kllz,r). The 
two open circles represent Tj^) and T2D calculated from 
Eq.^with d — 1 unit cell and d = film thickness, respec- 
tively. The thin solid curves are quadratic fits to ns(T) 
just below T2JJ. They approximate mean-field behavior. 
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FIG. 2: Xj^iT) = dflim/A^ and cri(50 kHz,T) vs. T measured 
at 50 kHz for underdoped YBCO film D (thick curves), 40 
unit cells thick. The center of the 2 K wide peak in o-i(50 
kHz,r) very nearly coincides with T2D (open circle). Inset 
shows an enlarged view of the transition region. 

AlJ(r), for T > Ti^. 

We expect to see 2D fluctuations in the 4-unit-ceU thick 
film just because it is so thin. Indeed, the top panel of 
Fig. 1 shows a rapid downturn in ns{T) oc AJ^(T), high- 
lighted by comparison with the quadratic fit (thin curve) . 
The drop occurs about midway between Tj*^, = 60 K and 
T2D = 70 K, so we cannot tell whether the effective thick- 
ness for fluctuations is one unit cell or the fllm thickness. 
Note that the downward curvature in X^{T) is essen- 
tially constant from 4 K to 60 K, meaning that TPFs are 
neghgible below T*jj = QOK. 

We interpret the drop in X^(T) as a vortex-pair un- 
binding transition partly because Tc very nearly coin- 
cides with T2D obtained from Eq^] (circles in Fig. 1), but 
also because the observed shift of cti and X'^ to higher 
temperature with increasing frequency is similar to 
behavior seen in a-MoGe fllms j21j . Transitions here are 
slightly broader than in a-MoGe fllms, presumably 
due to slight fllm inhomogeneity. If the downturn is not 
due to TPFs, then TPFs are even weaker than we con- 
clude here. 

In the 8 unit-ccU-thick fllm, the drop in ns occurs 
closer to T2D than T2*£3, suggesting that the effective 
thickness for TPFs is the film thickness, not a bilayer 
thickness. Using the quadratic fit as a proxy for mean- 
fleld behavior, we flnd that TPFs suppress Tc by less 
than 5 K in this film. The drop in 10 unit-cell-thick fllm 
is somewhat broader than in the other films. 

We now turn to the main focus of this paper, the role 
of TPFs in underdoped YBCO. Figure 2 shows Aj^(r) 
and (Ti{T) measured at 50 kHz for underdoped film D. 
From the fllms Tc of 34 K, we estimate its oxygen sto- 
ichiometry at Oqa, so it is strongly underdoped. The 



FIG. 3: A"^(r) for a single crystal of BSCCO The thin 
solid curve is a quadratic fit to A~^(r) between 60 K and 76 
K. The circle indicates T,*r, ~ 80 A". 



peak in cti is due to critical fluctuations; the peak value 
of (71 is more than 10^ larger than the fllms conductiv- 
ity above Tc- The width of the peak, about 2 K, is an 
upper limit on the inhomogeneity in Tc- Nothing in the 
data at T w T2JJ « 22 K indicates that fluctuations are 
strong. Instead of arcing downward, X'^{T) develops 
upward curvature between T^j^ and T2D, before flnally 
dropping at T2D ~ H-TK. Our view is that TPF effects 
are signiflcant only very close to Tc, and that they sup- 
press Tc only a few Kelvins in this strongly underdoped 
cuprate. This is our central flnding. 

Since we flnd TPFs to be weak in YBCO films, while 
very clean YBCO crystals show strong fluctuations, it is 
interesting to consider BSCCO, a cuprate in which inter- 
layer coupling is about 100 times weaker than in YBCO. 
Figure 3 shows measurements 113 of A-2(T) measured 
at 34.7 GHz on a high-quality optimally-doped BSCCO 
crystal. X"^{T) is linear in T at low T, and the low-T 
slope extrapolates to zero at about 125 K. Lee et al. as- 
signed Tc — from the peak in ai. The open circle 
in Fig. 3 marks Tj*^ w 80 K. The thin solid curve is a 
quadratic flt to data between 60 K and 76 K. Data drop 
below the fit for T > 88 K, which we interpret as the 
onset of strong phase fluctuations. We emphasize that 
at 88 K the properly normalized temperature, T/A^^(r) 
is 2.6 times larger than at Tj*^. We argue that TPFs 
must be weak at 80 K if superconductivity can survive a 
factor of 2.6 increase in normalized temperature. From 
the quadratic flt, or from a tangent to the data at 80 K, 
we estimate that TPF's suppress Tc by about 6 K. It 
seems that the fluctuations in very clean YBCO crystals 
are more interesting than simple phase fluctuations. 

Why are fluctuations so weak? It is known that TPFs 
weaken when T drops below an effective Debye temper- 
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ature, H, 0, but it turns out that this effect cannot 
account for our results. The effective Debye tempera- 
ture, Tq, is the temperature where the thermal frequency, 
ksT/h, equals the characteristic " R/L" superfluid relax- 
ation rate, ri(T) w piq / ^q\^{T). pjq is the normal-state 
resistivity. We estimate Tq « 45 K, 40 K, and 60 K, for d 
= 4, 8, and 10 unit cells, respectively, and Tq k, 25K for 
the underdoped film. Thus, in the optimally-doped films 
this quantum effect is certainly minor for T > Tj*^,. For 
the underdoped film, even though Tq lies between Tj*!? 
and T(7, quantum suppression of TPFs is by an algebraic 
factor, « + hn{T)/kBT], and would not be strong 
enough to obscure the dramatic TPF-induced drop inng. 

It seems that electron transport in cuprate supercon- 
ductors is somehow granular. Grains extend through the 
film thickness, even for films 40 unit cells thick, hence 
the effective thickness for phase fluctuations is the film 
thickness. Granular models have been proposed, e.g., 
ref. HI Hi El describe a model that accounts for many 
physical properties of cuprates, including the T- linear re- 
sistivity at optimal doping. 

Finall y, w e comment on fluctuation behavior in YBCO 
crystals [l^. In very clean YBCO samples, the 3D criti- 
cal region extends below Tj^,, whereas in the layered su- 
perconductor model the 3D critical region must lie above 
^2*D- F'o'^ example, Kamal et al. find that at temper- 
atures from - 85% to - 99.9% of Tc ~ 93.78 K, X-^{T) 
is fitted to within 1% by the 3D-XY exponent of 2/3: 
A-2(T)/A-2(0) « 1.26(1 - T/Tc)^/^. T*^ for this crys- 
tal is 0.96rc, where {T:^d) / i'^) ~ 0.14. Thus, the 
paradigm of cuprates as layered superconductors does not 
find good agreement with fluctuation behavior observed 
in the superfluid density. 

In conclusion, we have made high precision, low fre- 
quency measurements of \J^{T) and cri(T) in 4 to 40 
unit-cell thick YBCO films. Two-dimensional thermal 
phase fluctuation effects are observed, but the charac- 
teristic thickness is the film thickness, not the thickness 
of a Cu02 bilayer, and TPFs are therefore weaker than 
expected. Our main point is that even in strongly under- 
doped YBCO films, thermal phase fluctuations, as they 
are understood from simulations of Josephson-coupled 
grains, suppress Tc by only a few Kelvins. They cannot 
account for the reduction of Tc to zero with underdop- 
ing. It may well be that fluctuations of some kind sup- 
press the measured Tc well below its mean-field value, 
but that fluctuation is subtler than simple thermal phase 
fluctuations in a quasi-2D superconductor. 
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